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Abstract 



We propose a general method to obtain accurate estimates for some of the "low- 
energy constants" in the one-loop effective chiral Lagrangian by means of simulating 
^ ■ lattice QCD. In particular, the method is sensitive to those constants whose values are 

required to test the hypothesis of a massless up-quark. Initial tests performed in the 
quenched approximation confirm that good statistical precision can be achieved. As a 
byproduct we obtain an accurate estimate for the ratio of pseudoscalar decay constants, 
Fk/Fn, in the quenched approximation, which lies 10% below the experimental result. 
The quantities that serve to extract the low-energy constants also allow a test of the 
scaling behaviour of different discretizations of QCD and a search for the effects of 
dynamical quarks. 
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1 Introduction 



Chiral Perturbation Theory (ChPT) (T|,|2| plays an important role in the study of low- 
energy phenomena in QCD. As is well known, ChPT is based on simultaneous expan- 
sions of scattering amplitudes or hadronic matrix elements in powers of momentum and 
the quark masses m u , m& and m s , and the form of the effective chiral Lagrangian is 
entirely determined by chiral symmetry. Another feature is the emergence of a num- 
ber of coupling constants at every given order, which incorporate the non-perturbative 
dynamics of QCD. These couplings - sometimes called "low-energy constants" - are a 
priori undetermined and can be fixed through phenomenological constraints in conjunc- 
tion with additional assumptions. Currently the typical accuracy which is achieved in 
the determination of the low-energy constants is not very high (for a review, see Q). 

It has been noted some time ago that the mass parameters m u , m& and m s cannot 
be fixed unambiguously from symmetry considerations alone. The reason is that the 
effective chiral Lagrangian is invariant under a family of transformations of the quark 
masses, which can be absorbed into the low-energy constants Q|. In particular, this 
hidden symmetry of the chiral Lagrangian implies that one cannot distinguish between 
m u = and m n ~ 5MeV, whilst preserving all other phenomenological predictions. 
Thus, additional theoretical input beyond ChPT is required to decide whether m u = 
is indeed realized in nature. An obvious strategy is then to sharpen the current estimates 
for the low-energy constants and check whether they are compatible with the condition 
that m u = 0. 

Since m u = represents a simple solution to the strong CP problem, this question 
has continued to attract a lot of attention, despite several plausible arguments, each 
of which independently rules out m u = || ||. However, the problem has never been 
studied in a first principles approach starting from the QCD Lagrangian directly. 

In this paper we propose and test a method to determine a large set of low-energy 
constants with good accuracy using lattice simulations. Given its technical feasibility, 
such an approach eliminates the use of theoretical assumptions in the specification of 
the chiral Lagrangian. In the context of testing m u = the role of lattice calculations 
has recently been discussed by Cohen, Kaplan and Nelson |7}. We expand on their 
proposal by incorporating the information gained by studying the mass dependence 
of matrix elements in addition to that of the pseudoscalar masses. Furthermore, we 
present ready-to-use numerical procedures which increase the statistical precision and 
discuss the influence of lattice artefacts. The latter point is of great relevance because 
the effective chiral Lagrangian of Gasser and Leutwyler is not valid for non-zero lattice 
spacing. 

Our method is generally applicable in simulations of QCD with or without dynam- 
ical quarks. This initial study mainly serves to test its accuracy and limitations, and 
for that purpose it is sufficient to apply it to lattice data obtained in the quenched 
approximation. As a consequence we also consider comparisons of our lattice data with 
the quenched version of ChPT, thereby extracting some of the low-energy constants 
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relevant for the quenched theory. 

On the basis of our lattice data we conclude that the mass dependence of matrix 
elements can be determined with high precision in lattice simulations. Furthermore we 
show that low-energy constants in the chiral Lagrangian can be obtained with a typical, 
absolute statistical accuracy of ±0.05 in the continuum limit. Systematic uncertainties 
due to neglected higher orders are estimated to be ± 0.2. However, a systematic study 
of the influence of higher orders in the chiral expansions has so far not been performed, 
owing to limitations in the range of light quark masses that one is currently able to 
explore. Future calculations at smaller quark masses will be required in order to settle 
this point. 

The remainder of this paper is organized as follows. In Section ^ we recall the main 
concepts of ChPT. Our computational strategy is described in Section ||, and the numer- 
ical details are discussed in Section ^, focusing on the extrapolations to the continuum. 
Our results are presented in Section ||, and Section ^ contains some concluding remarks. 
Two appendices list some details about expressions in partially quenched ChPT and 
explain our choice of additional low-energy constants in the quenched theory. 

2 Chiral perturbation theory 

In this section we repeat the main features of ChPT and specify our notation. In addition 
to standard ChPT we also discuss its quenched and partially quenched versions. 

2.1 Standard ChPT and m u = 

The effective chiral Lagrangian is written as a low-energy expansion 

£ cS = £ ( $+£ ( $ + ... (2.1) 

(2) 

The leading contribution £ c{ { contains two coupling constants Fq and Bq. To lowest 
order Fq coincides with the pion decay constant F w . Throughout this paper we adopt 
a convention in which F n ~ 93MeV. At order p additional couplings ai,a2, ■ ■ ■ , U12 
appear in the effective chiral Lagrangian These couplings are not constrained by chiral 
symmetry. Furthermore they contain the divergences that arise if one goes beyond tree 
level and thus depend on the renormalization scale (and scheme). As mentioned in 
the introduction, experimental information at low energies is not sufficient to specify 
the values of the complete set of couplings a±, . . . , ctyi- For this reason one has to add 
further theoretical constraints, which usually involve certain assumptions, such as large- 
N c arguments. By applying such a procedure, the "standard" values of the low-energy 

2 We adopt a convention in which the aSs are related to the corresponding couplings Li of ref. Q 
through Qi = 8(47r) 2 Li. 
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constants in our convention for Nf = 3 flavours read Jl] [^,|| |llj 



a i = 


0.2 ±0.4 


a 6 = 


-0.5 ±0.4 


«2 = 


1.07 ±0.4 


CK7 = 


-0.5 ±0.25 


"3 = 


-4.4 ± 1.4 


"8 = 


0.76 ±0.4 


«4 = 


-0.76 ±0.6 


a 9 = 


7.8 ±0.8 


"5 = 


0.5 ±0.6 


UlO = 


-6.1 ±0.8. 



(2.2) 



Here the aj's have been renormalized at scale [i = 4:irF n , which will always be used in 
the remainder of this paper. The constants an and a\2 are of little phenomenological 
interest and are not included here. 

The question whether m n = has been discussed at length in the literature |@-@]. 
The usual starting point is the observation that simultaneous transformations of the 
quark masses 

= m n ± Xm^rris, m\ = + \m s m n , = m s + Xm n m^, (2-3) 

and coupling constants according to 

A , , (4vrF ) 2 A , , (4vrF ) 2 A N (4ttF ) 2 

a 6 =a 6 ±A— - — , a 7 =a 7 ±A— - — , a 8 = a 8 - A— - — , (2.4) 
4-Do 4±J 2i^o 

leaves the effective chiral Lagrangian invariant. Here, A denotes an arbitrary parameter, 
and for A / 0, m n = one can generate an effective up-quark mass such that all 
predictions by ChPT, in particular for ratios of quark masses, remain unchanged. In 
order to test whether m u = one has to determine the linear combination || (2a§ — 05), 
which, however, is not directly accessible from phenomenology. The value of 05 is 
obtained from the ratio of kaon to pion decay constants, Fx/F n , but as can only be 
inferred from the linear combination 

05 — 12a7 — 6ag, (2-5) 

which is related to the Gell-Mann-Okubo formula. Clearly this linear combination is 
invariant under the transformation of eq. (|2.4j ). As pointed out in a choice for a-j 
and as which is compatible with m u = is given by 

m u = 0: a 8 = -0.9±0.4, a 7 = 0.25 ± 0.25, (2.6) 

which is quite different from the corresponding numbers listed in eq. (pT2|). The task 



for lattice simulations is now to provide independent determinations of 05, as from 



linear combinations which are not invariant under eq. (2^), starting from the QCD 
Lagrangian. Provided that these estimates turn out to be sufficiently accurate, it should 
then be possible to test the hypothesis that m u = 0. 
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2.2 Partially quenched ChPT 

The role of lattice simulations for the determination of the c^'s has already been stressed 
in |Q,|l2|, and most recently in JDJ. In particular, it has been noted that simulations of 
"partially quenched" QCD, in which the sea and valence quarks can be chosen indepen- 
dently, may be useful as well as technically feasible. Thus, one is not forced to simulate 
at the physical values of the u, d and s quarks. Instead, the only requirement is that 
the pseudoscalar masses be small compared with the typical chiral scale of A x ~ 1 GeV. 
Hence the use of moderately light sea and valence quark masses and their independent 
variation may be sufficient to extract the low-energy constants. 

It has to be kept in mind, though, that the values of the a^s have to be determined 
separately for JVf = 2 and 3 flavours of dynamical quarks. This is a relevant point since 
there is currently not much experience with simulation algorithms for odd numbers of 
flavour. 

Partially quenched ChPT has been studied to 0(p 4 ) by a number of authors fl2l , 



14, 15 1 . Here we focus on the one-loop expressions for pseudoscalar masses and decay 
constants obtained by Sharpe [|l4| for Nt degenerate flavours of sea quarks. For the 
remainder of this paper we also take over some of the notation used in that reference. 
In particular, we denote the mass of the sea quark by ms, whereas the masses of the (in 
general non-degenerate) valence quarks are denoted by mi , m 2 ■ As in |L4|] we introduce 
the dimensionless parameters 



By setting A x = 4ttFq in eqs. (14) and (18) of ref. fl4|| , we obtain the generic formulae 
for the pseudoscalar mass and decay constant, i.e. 

yiiivss - yn)lnyn - V22(yss - 2/22)^12/22 



y 12 (4vrF ) 2 |l + ^- 
+yi 2 (2a 8 - ct 5 ) + yssN { (2ae - a 4 )| (2. 



2/22 - yn 



Fps N{ 

-pr- = 1 - -j- (Vis In yiS + 2/2S In 2/25) 

1 / yny22 - yssyu , yn . 

2iVf v 2/22 - yn 2/22 
1 Nf 

+yi2-^a 5 + yss^a A . (2.9) 

Here the constants on are to be evaluated at scale fi = A x = AttFq ~ 1170 MeV. These 
expressions will later be used to form quantities that allow for the extraction of the aj's 
using lattice data. 
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2.3 Quenched ChPT 



Quenched ChPT has originally been discussed in refs. [16, [17]. The complete chiral La- 
grangian to order p 4 in quenched QCD has been studied by Colangelo and Pallante [18]. 
Their results form the basis of our analysis. 

As is well known flavour singlets have to be treated differently in the quenched 
approximation: the decoupling of the rj' from the pseudoscalar octet by means of the 
anomaly is not realized in the quenched theory. Therefore, the quenched chiral La- 
grangian contains additional low-energy constants associated with flavour-singlet con- 
tributions. These include the singlet mass scale mo and the coupling constant a$, which 
multiplies the kinetic term of the singlet field in the quenched chiral Lagrangian. The 
mass scale uiq is usually expressed in terms of the parameter 5 defined by 



3(4vrF ( 



\2' 



(2.10) 



From various lattice studies (e.g. |19|-[22||) one expects 5 ~ 0.1. For a$ the available 
estimates have been summarized in [^3| as a§ ~ 0.6. 

Following ref. [^] we restrict ourselves to the case of degenerate quarks. The 
complete results at one loop for the pseudoscalar mass and decay constant readf] 

m 2 PS = y(4vrF ) 2 {l-(5-|a^)[l + lny]+((2^-^)-^y} (2.11) 

^ = 1 + vial (2.12) 

Here y is defined as 

and «g, ag denote the analogues of the low-energy constants as and a*, in the quenched 
theory. 

3 Strategy 

We now introduce the procedure to determine the low-energy constants from lattice 
data by studying the quark mass dependence of suitably defined ratios of pseudoscalar 
masses and matrix elements. In particular it is useful to investigate the mass dependence 
around some reference quark mass m re f. It is important to realize that this reference 



point does not have to coincide with a physical quark mass [25]. We only require 
that it should lie inside the range of simulated quark masses and within the domain of 
applicability of ChPT. 

3 In ref. a different combination of low-energy constants appears in the expression for m,p 3 , since 
the authors use an alternative convention for the counterterms [Q. The convention employed in this 
paper is consistent with that used in standard and partially quenched ChPT. 
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3.1 The ratios Ryi and Rp 

Let us consider the case of degenerate valence quarks, mi = m<i = m, either in the 
quenched approximation or in partially quenched QCD at a fixed value of the sea quark 
mass ms- If we introduce 

2B m 2B m rei 
then the ratios defined by 

Rm(x) = U^J/U^J (3 - 2) 

Rf(x) = F PS (y)/Fps{y rei ) (3.3) 

are universal functions of the parameter x, which measures the deviation from the 
reference point y re f. Here, Gps denotes the matrix element of the pseudoscalar density 
between a pseudoscalar state and the vacuum, and the arguments of Fps, Gps refer to 
the quark mass value at which the matrix elements have to be evaluated. Using the 
definition of the current quark mass in terms of i*ps, Gps and mps (see, for instance, 



eqs. (2.1) and (2.2) in ref. [25Q, eq. ( p^ ) can be rewritten as 



Ru(x) _ * ) (Jf-A = ,!&&f (3.4) 

Extracting the low-energy constants from the ratios Rx, X = M, F instead of the 
masses and matrix elements themselves has several advantages: 

• The ratios Rx can be computed with high statistical precision owing to the strong 
correlations between numerator and denominator; 

• The renormalization factors of the axial current and the pseudoscalar density 
drop out in Rm and i?F-Q The ratios can therefore be readily extrapolated to 
the continuum limit, so that all dependence on the lattice spacing is eliminated. 
Strictly speaking it is only in the continuum limit that one is justified to compare 
the predictions of ChPT with lattice data; 

• Since discretization errors are under good control in Rx the effects of dynamical 
quarks can be isolated unambiguously. 

The high level of statistical accuracy of the ratios is crucial in order not to compromise 
the precision in the continuum limit. Extracting the low-energy constants from the 



4 For O(a) improved Wilson fermions, there remains a small mass-dependent renormalization pro- 
portional to (6a — bp) Or — l)am le f and bA(x — l)ara rc f in the case of Rm and Rf, respectively. Our 
experience has shown |2q] that this can be safely neglected. In our calculations we set &a and bp equal 
to their one- loop perturbative values fed]. 
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x-dependence in the continuum limit in turn guarantees that these estimates will not 
be distorted by cutoff effects. 

The simple renormalization and high precision of the ratios may also be exploited 
to perform scaling tests for different fermionic discretizations. 



3.2 Expressions for Rm and Rp in ChPT 

Below we give a list with the expressions for Rm and i?F in quenched and partially 
quenched ChPT. For simplicity we restrict ourselves to the case of degenerate valence 
quarks. The case of non-degenerate valence quarks in partially quenched ChPT is 
discussed in detail in Appendix [A|. 

Note that we have not yet specified the reference point y re f . At this stage the precise 
definition of its numerical value is not needed, and we thus postpone its specification to 
Section ||, where we describe our numerical results. 

We begin by considering Rm and Rf in quenched ChPT. By inserting eq. ( [2,11 ) 
into eq. (3.4) we obtain 



1 - (5 - §a$2/ re f)(l + lny re f) + y re f [(2a§ ~ «e 



M 1 - (5 - §a<j>xy rcf )(l + ln(xy rcf )) + xy rci [(2a% - af) - ^] 

Provided that all masses and couplings are small, it is allowed to expand the denomi- 
nator, which gives 

Ryi(x) = 1 + <51n x — |a$2/ re f xlnx + (x — + ln?/ re f, 



Similarly we obtain 



-y„,(x - l)(2o§ - al). (3.6) 



_ i_t^M , (3. 7 ) 

1 + y r ef 2«5 

which, after expanding the denominator, becomes 

R*(x) = l + y re{ (x-l)^al (3.8) 
In partially quenched ChPT it is useful to always define the reference point y re f at 



2B m Te { 
(4vrF ) s 



ni! = m 2 = m s = m ref , y ref = //| _ r , 2 . (3.9) 



There are several possibilities to study the mass dependence of the ratios Rx- Let us 
first consider the case of equal sea and valence quark masses. The x-dependent part in 
Rx is then obtained by setting 

SS: mi = m-2 = ms = xm ie {, (3.10) 



7 



which will be labelled "SS" in the following. By taking the appropriate limits in eqs. (2.7) 
and (2.8) for the above choices of quark masses and inserting the resulting expressions 
into the definition of Rm, we obtain, after expanding the denominator: 



Rl!(x) 



1 



Nf 



y rci [x In x + (x - 1) In y ie{ ] 



-2/refO - 1) [(2«8 - as) + Nf(2a 6 - a 4 )] 
For Rp the corresponding result is 

Nf 



(3.11) 



Rf{x) 



1 - — y re{ [xlnx + (x - l)lny ref ] + y re{ (x - l)\ (a 5 + Nf 0*4) . (3.12) 



For mi, ni2 7^ ms the x-dependence can be mapped out using either the valence or the 
sea quarks. In the former case, which we label "VV" we define the x-dependent part 
through 



VV: mi = TU2 = xm rc f, ms = m re f 
instead of eq. ( |3.10| ), which leads to the expressions 



(3.13) 



< v (x) 



i? F vv (x) 



1 - ^^ref [(2x - 1) lnx + 2(x - 1) lny ref ] 



y Ic{ (x - 1) (2a 8 -a 5 ) + ^ 

(x + 1) In (i(x + 1)) + (x - 1) lny ref 



1 N f 



(3.14) 



(3.15) 



A comparison of the expressions for Rm and Rf for the "SS" and "VV" cases shows that 
they are sensitive to different linear combinations of low-energy constants, depending on 
whether the x-dependence is defined using eq. (|3.10[) or eq. ( p.l3|) . In particular, we find 
that it is possible to extract directly from R M ^(x) the linear combination (2ag — «s), 
which is relevant to the question of whether m u = 0. 

There are several other possibilities to define the dependence of Rx on the quark 
masses, also allowing for non-degenerate valence quarks. Details are listed in Ap- 
pendix |A|. 



3.3 Extracting the low-energy constants 

We now describe a method to determine the low-energy constants from the ratios Rx 
in the quenched and unquenched cases. To this end we choose two values of mass 
parameters, x±, X2, and introduce the quantity 

AR x (x 1 ,x 2 ) = Rx(xi) - Rx(x 2 ), xi<x 2 , X = M, F. (3.16) 
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By inserting the expressions for the ratios Rx we can easily solve for the low-energy 
constants. For instance, from eqs. ( |3.8j ) and ( |3.16| ) we find 



i AR^(xi,x 2 ) 
'(xi - x 2 )y iei ' 



a* = 2 ~ FK yC . ( 3 - 17 ) 



Similar combinations of the a^'s can be obtained in partially quenched QCD. The explicit 
expressions are given in Appendix [A], and one can easily convince oneself that they serve 
to determine a^, 05, uq and as- 

The quantity AR^ is a special case, since it also depends on the low-energy con- 
stants 5 and q$, which only occur in the quenched theory. However, by inserting the 
estimates for S and a<j> quoted in the literature we can solve for (2«g — a§), i.e. 

(2o$-o£) = {y re{ (x 1 - X2)}' 1 x ^5ln(x 1 /x 2 ) - ARl l (x 1 ,x 2 ) 



xilnxi -x 2 lnx 2 + (xi - x 2 ) (\ + m y re f) }• (3.18) 



The expressions for Rx themselves can also be used in order to extract the low-energy 
constants from least-^ 2 fits over a suitably chosen interval in x. The differences ARx, 
however, have the advantage that some of the lattice artefacts may cancel. Thus, 
instead of first extrapolating Rx to the continuum limit and then forming the differences 
ARx, one may compute Ai?x(xi,X2) at non-zero lattice spacing and then perform the 
continuum extrapolation. Obviously the results must be independent of the ordering 
of the two procedures, which offers an additional check against the influence of lattice 
artefacts. By contrast, if a fitting procedure is employed at non-zero lattice spacing 
then it is a priori not clear whether the continuum expressions for the ratios Rx are 
valid. 

Our method to extract the low-energy constants is only viable if there is sufficient 
overlap between the region of validity of ChPT and the range of quark masses accessible 
in current lattice simulations. On present computers it is not possible to simulate very 
light quarks without suffering from large finite-volume effects. Furthermore, the fermion 
action and lattice spacings employed in this work do not allow the use of quark masses 
which are below approximately half of the strange quark mass. The reason is the 



occurrence of "exceptional configurations" [ 27 , 28 1 , which correspond to unphysical zero 



modes of the lattice Dirac operator. Therefore, since one is restricted to a range of 
relatively large quark masses, one must check the results against the influence of higher 
orders in the chiral expansion: if large, these would modify the numerical estimates for 
the low-energy constants considerably. This will be discussed in more detail in Section ||. 

4 Continuum limit of Rx and ARx 

The ratios Rx have been computed using the same quenched configurations as in our 



earlier papers [25, 29j. Results for hadron and current quark masses, as well as for Fp§ 



and Fps/Gps are listed in Table 1 of ref. [25]. Details of our numerical procedures are 
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Rm(x) 



R F (x) 



1.004 
1 

0.996 
1.004 
1 

0.996 
1 

0.996 
0.992 
0.988 



1.02 



x=1.10 



x = 0.90 




x = 0.75 



_L 



i - 



0.01 0.02 
(a/r ) a 



0.03 



1.01 
0.99 



0.98 



0.96 



0.95 



x=1.10 



x = 0.90 





— f— — ' i - 




x = 0.75 
1 , , , , 



0.01 0.02 0.03 

(a/r )2 



Figure 1: Continuum extrapolations of Rm(x) and Rf(x) for selected values of the mass parameter x. 



described in Appendix A of the same paper. Here we only mention that non-perturbative 
O(a) improvement |27,3C,3l| has been employed, and we will therefore assume that the 
remaining discretization errors are of order a 2 . 

In contrast to our earlier papers the statistical errors in this work have been es- 
timated using a bootstrap procedure [32|, with 200 bootstrap samples generated from 
the sets of gauge configurations. This allows us to keep a constant number of bootstrap 
samples at every value of the bare coupling, regardless of the number of configurations. 
By performing the continuum extrapolation of Rx(x) for every bootstrap sample, our 
error procedure thus preserves the correlations in the mass parameter x. Through- 
out this paper we quote the symmetrized error from the central 68% of the bootstrap 
distribution. 

The value of the reference quark mass m re f is defined through the condition 



(mp S r ) 



3, 



(4.1) 



a choice that has also been considered in refs. [25, 33 1 . For Fq = F n = 93.3 MeV the nu- 
merical value of y re f is thus y re f = 0.3398 .... Lattice data for the hadronic radius ro [34] 
have been taken from ref. |35|]. For ro = 0.5 fm the definition of eq. ( |4.lD corresponds 
to a pseudoscalar meson whose squared mass is roughly twice as large as the kaon mass 
squared, and therefore m re f ~ m s (with "s" standing for "strange"). The results for Fp$ 
and -Fps/Gps obtained for several values of the bare coupling have been interpolated 
in the current quark mass m to common values of x = m/m re f. To this end the two 
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X 



-0.025 
-0.03 



K -0.035 



-0.04 


* -0.005 



K -0.01 
< 



-0.015 



i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r 



; Xj = 0.75, x 2 = 0.95 

J I I I i i I i i i i I i i i i I i i i u 



— i 1 1 1 1 1 1 1 1 1 r 



x^O.75, x 2 = 0.95 



J I I I I I I I I I I I I I I I I I I I L 



0.01 0.02 
(a/r )» 



0.03 



Figure 2: Continuum extrapolations of A_Rm and ARp. 



neighbouring points which straddle the value of x were used in a linear interpolation. 
If x was slightly outside the range of simulated quark masses, a linear extrapolation was 
performed using the three nearest points. The stability of this procedure was checked 
by varying the input masses and/or performing quadratic interpolations/extrapolations. 
Our sets of simulated quark masses cover the range 0.75 < x < 1.4, and we have chosen 
increments of 0.05 to map out the mass dependence of Rx{x). 

In Fig. |l] we plot the ratios Rm(x) and Rp(x) against (a/ro) 2 for three different 
values of x. The plots show that lattice artefacts are very small in general and are 
consistent with leading cutoff effects of order a 2 . As a safeguard against higher-order 
lattice artefacts, we have excluded the points computed for our coarsest lattice (a « 
0.1 fm) from the extrapolation. Results obtained by performing the extrapolations for all 
four values of the lattice spacing were entirely consistent, albeit with a smaller statistical 
error. 

At non-zero lattice spacing the statistical precision is better than 0.3% and 0.2% for 
Rm and Rp, respectively. In the continuum limit these figures are only slightly larger, 
namely 0.4% for Rm and 0.3% for Rp. This demonstrates the high level of precision that 
can be achieved in the continuum limit. Furthermore it is clear that the extrapolation 
is well controlled. 

As mentioned in Section ||| the results for ARx in the continuum limit can be 
obtained either directly from the continuum values of Rx or from a continuum extrap- 
olation of ARx computed at non-zero lattice spacing. The latter extrapolations are 
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shown in Fig. |2] for a particular choice of x\ and X2- As in the case of the ratios Rx 
themselves the continuum limit is easy to control. 



5 Results 

We can now compare our results for the ratios Rx to the expressions predicted by 
ChPT. Since the numerical data have been obtained in the quenched approximation, 
we will focus on the determination of 05 and (2a§ — a§). 

5.1 Low-energy constants in quenched ChPT 



The determination of from eq. ( [3.17 ) is straightforward, since it only depends on the 



mass parameters x\ and x-i- However, in order to compute (2ag — a§) from eq. ( |3,18|) 
one must make a suitable choice of 6 and a$. Here we are going to consider two cases, 
namely 

Ql: 5 = 0.12 ±0.02, a$ = (5.1) 
Q2: S = 0.05 ±0.02, a$ = 0.5. (5.2) 

The reasoning which led us to consider these two distinct sets of parameters is described 
in Appendix [b]. 

Our estimates for ot§ and (2«g — ot^) have been obtained from ARp(x±, X2) and 
A.Rm(x\, X2) for x\ = 0.75 and X2 = 0.95. This choice was motivated by the desire to go 
to the smallest possible quark masses, whilst maintaining a reasonably large x-interval 
in order to check the stability of the results against variations in the parameters x\ 
and X2- Following this procedure, we obtained the following estimates for the low- 
energy constants in quenched ChPT 

al = 0.99 ±0.06 (5.3) 

q q / 0.35 ± 0.05 ± 0.07; Ql 

(2a «-^ } = \ 0.02 ± 0.05 ± 0.07; Q2 " (5 " 4) 

Here, the first error is statistical, while the second (where quoted) is due to the variation 
of ±0.02 in the value of 6 for both parameter sets Ql and Q2. These results can now be 
inserted into the expressions for the ratios Rm and Rf, eqs. (|3.6|) and (|3.8[) , respectively, 
and the resulting curves are compared with the data in Fig. ||] (a)-(c). 

The linearity of Rf predicted by one-loop quenched ChPT is well reproduced by the 
numerical data, resulting in very stable estimates for a^. Although the result in eq. (|5.3[ ) 
has been extracted from a fairly narrow cc-interval, it provides a good representation of 
the data over the entire range of quark masses considered here. The determination of 
(2ag — CK5) is also quite stable for both sets of values of 5, a<j>, i.e. Ql and Q2. In both 
cases a good description of the numerical data is achieved for x < 1.2, with (2«g — a^) 
extracted for Ql tracing the data quite well even at the largest x- values. 
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Figure 3: (a): Data for the ratio Rf compared with the curve obtained by inserting the low-en ergy 
constant of eq. (5.5); (b) and (c): the same for the ratio Rm using the low-energy constants of eq. (5.4) 
for parameter sets Ql and Q2, respectively. Dotted lines indicate the variation due to the statistical 
uncertainty in the extraction of the low-energy constants. 



Estimates for and (2a§ — a?) obtained from the continuum extrapolations of 
ARp and Ai?M shown in Fig. ^ are entirely consistent with eqs. ( |5.3D and (|5.4j) . The 
same is true if the low-energy constants are extracted directly from fits to Rx(x) for 
x in the interval 0.75 < x < 0.95; the results are numerically almost identical to those 
obtained using the differences ARx- 

5.2 Effects of higher orders in the quark mass 

Although the results presented in the last subsection suggest that the one-loop formulae 
for quenched ChPT offer a good description of the x-dependence of the ratios Rx, 
further work is needed to corroborate these findings by extending the range of quark 
masses under study to smaller values. For instance, the downward curvature observed 
in the prediction for Rm when x < 0.75 remains to be verified. 

Furthermore, data at smaller quark masses will be required to systematically assess 
the influence of higher-order terms in the chiral expansion on the determination of 
the low-energy constants. Such terms manifest themselves in additional contributions 
proportional to x 2 in the formulae for Rx- Since the range 0.75 < x < 0.95 corresponds 
to pseudoscalar meson masses between 590 and 670 MeV, it cannot be excluded that 
higher-order terms have a sizeable impact on the extraction of the a^s. 

Even without results at smaller masses there are several possibilities to study the 
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effects of neglecting higher orders in ChPT on our estimates for 05 and (2o§ — ctg). We 
stress, however, that none of the methods described below can replace the systematic 
investigation of smaller quark masses. 

An obvious way to proceed is to add a quadratic term to the expressions for R F 
and R%. For instance, the modified expression for R F reads 

R F (x) — » 1 + yref(x - l)\al + y 2 ef (x 2 - 1) d F , (5.5) 

with a similar quadratic term proportional to g?m m the case of R^. One can now 
perform least-x 2 fits over the entire range 0.75 < x < 1.4, thereby extracting the low- 
energy constants as well as dp, dyi- Because of the linearity of the data for Rp the only 
effect of including x 2 -terms in the determination of 05 is to increase its statistical error 



to ± 0.19. The estimates for (203—05) are more sensitive: compared with eqs. (5.4) their 
central values increase by 0.11 and 0.22 for the parameter sets Ql and Q2, respectively, 
while the statistical error is increased to ±0.12. The variation in the central values, or 
the larger statistical error in the case of Og, may serve as an estimate of the uncertainty 
induced by neglecting higher orders. 

An alternative, albeit naive, estimate of the effect in question is obtained by assum- 
ing that the chiral expansion converges like a geometric series. This implies that one 
expects the quadratic contributions to Rx to be roughly as large as the square of the lin- 
ear ones. Here we consider R F as the generic case, since it does not involve logarithmic 
terms. Its linear contribution amounts to ~ 16% at x = 1, so that the quadratic term 
is estimated as 0.025 a; 2 . If we generalize this estimate, then the systematic uncertainty 
in ARx(x±,X2) due to neglecting quadratic terms is given by 

syst. error in ARx{x\, X2) ~ ± 0.025 (x 2 — x\). (5-6) 

Through eqs. ( |3.17 ) and ( 3.1 8| ) this is easily translated into systematic errors^] on 05 

4 - <4) 

syst. error in Og ~ ± 0.25 (5-7) 
syst. error in (2o§ — 05) rs ±0.13, (5-8) 



and (2og — 05), as 



which is of the same order of magnitude as the previous estimate. 

Finally one can compare the expanded and unexpanded expressions for Rm and Rp , 
which differ at order x 2 (cf. eqs. (3.5)-([T8|)). By extracting the low-energy constants 



from least-x 2 fits to the ratios in eqs. ( |3.5|) and ( |3.7|) we obtain yet another set of results. 
For 05 the central value is larger by 20%, whereas the result for (2o§ — Og) remains 
essentially unchanged. 

In order to present a synthesis of the various efforts to estimate the uncertainty 
due to neglecting higher orders, we note that the typical size of this systematic error 
amounts to ±0.2 for both 05 and (2o§ — 05). 

5 Instead of the constant 0.025 in eq. (^^) the reader may supply an alternative value, depending on 
whether our estimate is deemed too optimistic or pessimistic. 
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5.3 Application: the ratio F^/F^ 

The result for a§ extracted from Ftp can be used to compute the ratio of the decay con- 
stants of the kaon and pion, F^/F^. In fact, one usually employs the reverse procedure 
by using the experimental result for F-^/F^ to extract the phenomenological value of 

If we assume that contributions proportional to differences in the valence quark 
masses can be neglected, we can simply use the definition of Rp(x) to compute F^/F^: 

IT = pq£ K | = 1 + y™t( XK ~ x ^)\ a b- ( 5 - 9 ) 

^■n -ftp \X-ir } 

Here the dimensionless mass parameters xk and x n are related to the kaon and pion 
masses by 



2 

m 



XK Vvei 



{^F Q f 

2 

mi 



§^ = 0.1782 (5.10) 



s*yref = (A " 2 =0-0142, (5.11) 
(4vrF )^ 



where, as in ref. [25], we have used twk = 495 MeV and m n = 139.6 MeV. If the estimate 
for «g from eq. (5.3) is inserted we find 

Fk 

= 1.081 ±0.005 ±0.017. (5.12) 

Fir 

Here the first error is statistical and the second is the estimated uncertainty in ot^, 
due to neglecting quadratic terms. The above result is significantly smaller than the 
experimental value of Fk/F w = 1.22 [36]. 

A formula for Fyi/F n in one-loop quenched ChPT, which also accounts for differ- 
ences in the valence quark masses, can be derived from eqs. (18) and (20) of ref. [37]. 
The low-energy constant that appears in the one-loop counterterm is the same as in 
the quenched degenerate case [18, 24 . In our notation the full one-loop expression for 
Fk/F-k reads 



^ = 1 + y ref (x K - x n ) - - - | (s - ^ x K y vei 



35 x K - a$ ^ re fX 7r (2x K - x n ) ^ / 2x K -x„ y (5 13) 

Note that one recovers eq. (|5.9| ) when 5 = a$ = 0. For the two sets of parameters, Ql 
and Q2, the results for F^/F^ are evaluated as 

Fk _ j 1.125 ± 0.005 ± 0.016 ± 0.007; Ql 

F n ~~ i 1.110 ±0.005 ±0.017 ±0.007; Q2 ' ^ ' ' 
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where the additional third error is due to the variation of ±0.02 in the input value 
for 5. While contributions proportional to differences in the quark masses enhance the 
result for Fk/F w compared with eq. ( |5.12 ), the value is still smaller than experiment by 



about 10%. Deviations of this order of magnitude are typical of quantities computed 
in the quenched approximation. This has previously been inferred from calculations 
of the hadron spectrum [jjj, and the results presented here firmly establish these 
findings for matrix elements of local currents as well. The fact that Fyi/F n is typically 
underestimated in quenched calculations has been observed before (3^,^]. Note, how- 
ever, that our estimates have much smaller uncertainties than those quoted in the other 
references. 

Finally we remark that the enhancement in the estimate of F^/F n due to differences 
in the quark masses, demonstrates that these effects can be quite significant if the quark 
mass difference is as large as that between the physical light and strange quarks. By 
contrast, estimates of these effects based on masses in the region of that of the strange 
quark tend to be much smaller p5| , |40[ ] . 

5.4 Partially quenched ChPT 

In addition to our analysis of the quenched version of ChPT we can also tentatively use 
the expressions for Rm and Rp which are derived in the partially quenched theory. This 
essentially serves two purposes. On the one hand we can investigate to what extent 
the ratios Rx in quenched QCD are described by the formulae of partially quenched 
ChPT. If a severe mismatch is encountered it will signal that sea-quark effects are 
not accounted for. On the other hand, by extracting the c^'s using the expressions of 
partially quenched ChPT but quenched numerical data, we can test how much of the 
relevant physical information for the low-energy constants in full QCD is encoded in the 
mass dependence obtained in the quenched approximation. The subsequent calculation 
of quantities such as the ratio Fyi/F n and its comparison with experiment provides a 
quantitative criterion for this task. In this spirit we have investigated the cases labelled 
"SS" and "VV" discussed in Section |. 

If we employ the expression for Rp S , eq. ( 3.12 ), we find that the linearity of the 
numerical data is incompatible with the additional logarithmic terms contained in the 
formula. As a consequence an acceptable representation of the data for Rp cannot 



be achieved on the basis of eq. ( |3.12| ). We conclude that the mass dependence of the 
quenched decay constant is incompatible with the chiral behaviour predicted in full 
QCD! 

By contrast, the expression for Rf^ fits the numerical data very well over almost the 
entire range in x. However, the disagreement between the data and i?| s indicates that 
the "SS" formulae cannot be applied to the quenched data in any reasonable manner. 
Therefore, we refrain from quoting an estimate for (2«s — 05) + 3(2a6 — 0:4), even though 
a good modelling of the data for Ryi can be achieved. 

On the other hand, the expressions for both ratios obtained for the "VV" case do 
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Figure 4: The "VV" case in partially quenched ChPT (a): data for the ratio Rf compared with 
the curve obtained by inserting the low-energy constant of eq. (5.15); (b): the same for the ratio _Rm 
using the low-energy constants of eq. (5.16). Dotted lines indicate the variation due to the statistical 
uncertainty in the extraction of the low-energy constants. 



indeed lead to an acceptable representation of the data, at least for x < 1. Quoting 
only statistical errors, the results for the low-energy constants read 

a 5 = 0.75 ±0.06 (5.15) 
(2a 8 -a 5 ) = 0.15 ±0.05. (5.16) 

The curves corresponding to these values are shown in Fig. || (a) and (b). While for 
x < 1 the curves fit the data quite well, they deviate much more at larger x than the 
corresponding quenched ones. Compared with the quenched case, the results for the at 
are roughly in the same ball park. 

Using the value for 05 extracted from i?p V we can now compute the ratio F^/F^. 
Since the low-energy constants in the partially quenched and full theories are the same, 
we have used the expression derived in full QCD H. In our notation it reads 



F k 

+ j( 32; 7r - 2x k) in(x n y rc{ ) - x K ln(x K y: 
-i(4x K -x.)ln 



1 + 2/ref(^K - X n )\a^ 



VK-x^lnf 43 * X " \\. (5.17) 
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After inserting the numerical values for xr, x w , y re f and «5 one finds 

= 1.228 ± 0.005 (stat) ± 0.016. (5.18) 

Here the uncertainty of ± 0.016 is due to neglecting higher orders. This result is actu- 
ally compatible with the experimental value, which is not entirely surprising, since our 
estimate of 05 = 0.75(6) agrees quite well with its phenomenological value of 0.5 ± 0.6 
quoted in eq. ( |2.2| ). These findings suggest that the valence quark mass dependence 
of the numerically determined ratios Rx in the region of the strange quark mass is 
only weakly distorted by neglecting dynamical quarks. The relevant quark mass effects 



which account for the difference between the quenched result eq. ( 5.14 ) and the exper- 
imental value Fk/Ft,- = 1.22, appear to be due to the very light up and down quarks 
in full QCD. Their non-linear effects can be efficiently accounted for by the formulae of 
partially quenched ChPT. 



6 Summary and conclusions 

In this paper we have proposed a general method to extract the low-energy constants in 
effective chiral Lagrangians by studying the mass dependence of suitably defined ratios 
of matrix elements and matching it to ChPT. The ratios are typically obtained with high 
statistical precision and can be extrapolated to the continuum limit in a straightforward 
manner. Thus, the method offers a conceptually clean way to determine the low-energy 
constants, since it is strictly speaking only in the continuum limit that a comparison of 
lattice data to ChPT can be performed. 

The dimensionless ratios Rx are also interesting in their own right: they can be 
used to study the scaling behaviour of different fermionic discretizations without having 
to compute renormalization factors. Furthermore, the effects of dynamical quarks can 
be isolated unambiguously, since lattice artefacts are under good control. 

In this initial study we have presented results for the quenched approximation. 
The typical value of the absolute statistical error of ± 0.05 in the low-energy constants 
confirms that the achievable precision is indeed quite high compared with "conventional" 
phenomenological determinations. The main limitation at present is that results at 
smaller quark masses are not available. It is therefore not possible to systematically 
investigate whether there is sufficient overlap between the range of simulated quark 
masses and the domain of validity of ChPT. As a consequence the present estimated 
uncertainty of ± 0.2 in the low-energy constants due to higher-order terms are fairly 
coarse. It is expected, though, that the simulation of smaller quark masses will be 
feasible within the formulation of QCD with a "twisted" mass term (tmQCD) [^T|. In 
this construction the mass parameter protects against unphysical zero modes of the 
Dirac operator, which alleviates the problem of exceptional configurations encountered 
in simulations with Wilson fermions. Initial studies have shown promising results, and 
it is planned to exploit this method further in the present context. 
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Returning to the problem of examining the question of whether m u = 0, it is useful 
to recall that the problem can be solved through a reliable determination of ag. In the 
quenched approximation we obtain 

q _ f 0.67 ±0.04 (stat) ±0.04; Ql 
° 8 ~ \ 0.50 ± 0.04 (stat) ±0.04; Q2 ' ( } 

where the second error is due to the variation in S, and the additional uncertainty arising 
from neglecting higher orders is estimated to be ± 0.2. There is a priori no reason why 
the low-energy constant Qg should be numerically close to its counterpart in the full 
theory, and hence it would surely be premature to give a full assessment of the problem 
on the basis of our quenched results. Nevertheless, it is interesting to note that our 
results for a§ are quite close to the "standard" value for a$, eq. fl2.2|), which supports a 
non-vanishing up-quark mass. By contrast, given our data for the ratios Rx it would be 
very difficult to accommodate a large negative value, which is required for the up-quark 
to be massless (see eq. (|2.6|) ). 

Finally we wish to add a few general comments on the philosophy of our approach. 
Provided that the quark masses used in simulations are small enough for the one-loop 
expressions of ChPT to apply, our procedures altogether avoid chiral extrapolations of 
lattice results, which are known to be quite difficult to control. The example of the 
computation of F-^/F^ shows that the problem can be separated in two parts. The 
information which involves quark masses near the chiral limit can be extracted safely 
in ChPT, whereas the mass dependence for masses in the region somewhat below the 
strange quark mass can be adequately studied in lattice QCD. The latter part is the ad- 
ditional theoretical input needed to determine some of the low-energy constants, which 
are not accessible from chiral symmetry considerations alone. In short, our method 
amounts to exploiting the complementary character of ChPT and lattice QCD. The 
applicability of ChPT itself can be tested once quark masses somewhat smaller than 
those in our present work become accessible. 

The next step is the application of our method to the case of dynamical quarks. 
Here the relevant formulae that will be required are listed in Appendix |Aj Simula- 
tions using either the O(a) improved Wilson action for Nf = 2 flavours of dynamical 
quarks or other improvement schemes [ 43 1 have already been performed [44|-47|. 



Whereas an analysis of those results will be able to treat the two-flavour case, the ex- 
tension to the physically most interesting case of iVf = 3 flavours will require a significant 
amount of additional simulations. It is also worth investigating applications to the case 



of non-degenerate sea quarks and to try to determine the low-energy constant aj [13]. 
However, the main message of this paper is the following. In order to settle the question 
of whether m u = one does not even require the same level of accuracy as that of the 
quenched results presented here. A satisfactory analysis of the problem in the case of 
three degenerate flavours is therefore quite a realistic prospect. 
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Pu(x;y Te {) A M (a) 



ss 


x In x + [x — 1) In y re f 


(2a 8 - a 5 ) + iVf (2a 6 - a 4 ) 


vv 


(2x - 1) In x + 2(x - 1) In y rcf 


(2a 8 - as) + i 


SS2 


-(a? - l)lny ref 


iV f (2a 6 - 04) - ^ 


VS1 


x In x + (x — 1) In y re f 


^(2a 8 - as) 


VS2 





i(2a 8 - as) + Nf(2a 6 - a 4 ) 



Table 1: The functions p and A for the ratio Rm(x). 



A Partially quenched ChPT — the general case 

In this appendix we list the expressions for Rm and Rp, as well as for ARm and ARp, 
in partially quenched QCD. In particular we discuss all possibilities to define the in- 
dependence of Rx, allowing also for non-degenerate valence quarks. 

We start by considering the expressions in eqs. ( ]2.7| )^( |2~9| ), At the reference point 
we require all quark masses to coincide with m re f (cf. eq. Q3.9p ). Two cases, labelled 
"SS" and "VV" have already been discussed in Section |3.2| , In addition to "VV" , one 
can also study the x-dependence by varying the sea quark mass ms for fixed, degenerate 
valence quarks: 

SS2: mi = = Tn Te f, ms = xm re f. (A.l) 

For non-degenerate valence quarks one may define 

VS1: mi = xm re {, m<i = ms = m re f (A. 2) 

VS2: nil = wirefi m i = m s = xm re f. (A. 3) 

In order to list the expressions for Rm and Rp for all cases SS, VV, . . . , VS2 it is con- 
venient to introduce the general parameterization 

R M (x) = 1 - ^pm(x; j/ref) - y rcf (x - l)A M (a) (A.4) 
R F (x) = 1 - Nfy rci p F (x;y Tc{ ) - y Tc{ (x - l)A F (a). (A.5) 

Here, pyi and p-p are functions of x and y re f, and Am, Af denote linear combinations of 
the low-energy constants. The expressions for p and A are shown in Tables |l] and ^, 
respectively. 

Using the functions px and Ax, X = M, F, it is now quite easy to solve for 
particular linear combinations of low-energy constants by considering the differences 
ARx(xi,X2) introduced in eq. ( p. 16 ). 
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PF{x;y re i) A F (a) 

SS \ (xlnx + (x — 1) lny re f) — |(as + iVfC^) 

VV l(-±iln^i + ^ilny ref ) -±a 5 
SS2 I (-±iln^±i + ^ilny rcf ) -fa 4 
VS1 i ^ln^ + ^lny ref + ^lnz) -|(a 5 + ^) 

VS2 \ ^ln^ + ^llny vei + ^xlnx^ -jag - f« 4 + 4^ 
Table 2: The functions p and A for the ratio i?p(:r). 

With these definitions the linear combinations of low-energy constants denoted by 
Xx(a), X = M, F are simply given by 

AR M (X1,X 2 ) + %f (pM(>l,y re f) - PM(>2,yrcf)) s 

A M a = ^— p r A.6 

y rof [xi - x 2 ) 

AR F {x 1 ,x 2 ) + y Ie fN i (pF(xi,y Ie f)-pF(x2,y ie f)) , . 

A F (a) = r . A.7 

y re f (xi - x 2 ) 

By evaluating the right-hand side, using numerical data for ARx, one can easily solve 
for the desired combination of a^s. 

B Choices for 5 and a$ 

In this appendix we motivate our choices for the parameters 8 and a$ specified in 
eqs. ( |5.lD and (fT^). Here we will closely follow the analysis presented by the CP-PACS 
Collaboration l2ltl, i.e. we consider the ratio 



ml s ( mi ,m 2 ) 2m, m^rm.ma) 2m 2 

s — ' — 9 — 7 r ^ ■ — 9 — 7 \> k**'*-) 

mi + m 2 mp S (mi,mij nil + m 2 mp S (m 2 ,m 2 ) 

where the arguments in irip S have been included in order to distinguish between degen- 
erate and non-degenerate mesons. By inserting the expressions for mp S in quenched 
ChPT for degenerate and non-degenerate quarks (cf. eq. (9) in |l6|]) one obtains after 
expanding the denominator^ 

1 , x fo 2/H + 2/22 , yn \ . / , . x . 2/22 , yn \ m , 
s = 1 + 5 2 In — + — -(yn + y 22 ) + 2yn In — . (B.2) 

V 2/11-2/22 2/22/ o V 2/11-2/22 2/22/ 

In ref. [21] the quantity t was defined as 

f = 2 _ 2/11 +2/22 lQ y_n_ (B3) 
yn - y22 y22 



'Note that the one-loop counterterms drop out in the expanded version of s. 
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-0.6 



-0.4 -0.2 
t 







Figure 5: The "wedges" defined through eq. (B.2) for different vaiues of S and a<s> as expiained in the 
text. 



For a$ = eq. ( [B.2D reduces to a simple, linear relation between s and t: 

s = l + 6-t, (B.4) 

with the slope given by 5. By plotting the numerically determined values for the ratio s 
against t for a<j> = 0, CP-PACS concluded that all their data were enclosed within the 
"wedge" defined by 6 = 0.08-0.12. This wedge is represented in Fig. || by the solid lines. 

So far the parameter a<j> has not been considered in this kind of analysis^]. We are 
now going to argue that a value of a<j> « 0.5, as suggested by Sharpe in [^J, is only 
compatible with the CP-PACS data enclosed by the solid wedge, if 5 is chosen in the 
range 0.03-0.07. To this end we identify mi = m s , so that yn = 0.343. The mass 
ratio yn/2/22 is then varied between 1 and 24.4, where the latter number is the value of 



2m s /(m u + TTid ) computed in standard ChPT 48]. For 6 = 0.08-0.12, a$ = one thus 
recovers the wedge denoted by the solid line in Fig. [j| which encloses the CP-PACS 
data. 

If one combines 6 = 0.08-0.12 with the independent estimate of a$ = 0.5 from |23[ ] 
one obtains instead the wedge defined by the dotted lines, which is clearly incompatible 
with CP-PACS's numerical data. Agreement can only be restored if 5 is lowered. The 
dashed curves, which correspond to 5 = 0.03-0.07, a$ = 0.5, are again consistent with 
the wedge defined by the solid lines. These observations lead us to consider the two sets 



of parameters as specified in eqs. (5.1) and (157 



7 Terms proportional to q$ were also neglected in ref. [E2|, where 8 pa 0.08 was quoted. 
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